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We introduce a circuit quantum electrodynamical setup for a quantum single photon transistor. 
In our approach single photons propagate in two open transmission lines that are coupled via two 
interacting transmon qubits. The interaction is such that photons are not exchanged between the two 
transmission lines but a photon in one line can completely block respectively enable the propagation 
of photons in the other line. High on-off ratios can be achieved for feasible experimental parameters. 
Our approach is inherently scalable as all photon pulses can have the same pulse shape and carrier 
frequency such that output signals of one transistor can be input signals for a consecutive transistor. 

PACS numbers: 42.50.Ex,85.25.Cp,42.25.Fx,42.50.Nn 



Photons are the most suitable carrier for transmitting 
information over long distances as they are largely im- 
mune to environmental perturbations, and can propa- 
gate with very low loss and long-lived coherence in a 
wide range of media pQ. The use of photons in infor- 
mation processing however still suffers from the inability 
to realize controlled, strong interactions between indi- 
vidual photons. To make photons a more versatile in- 
formation carrier, it is therefore of great importance to 
conceive means of making photonic signals interact with 
each other [2]. In vacuum, direct photon-photon inter- 
actions are absent. Nonetheless, optical signals can in- 
fluence each other in nonlinear media. Yet, the quantum 
regime with interactions between individual photons only 
becomes accessible for devices where optical nonlineari- 
ties exceed incoherent and dissipative processes. Such 
devices therefore require a strong coupling of the pho- 
tons to the material that mediates the effective photon- 
photon interactions. Since the coupling of light to matter 
can be enhanced if light fields are confined to small vol- 
umes in space, cavities and one-dimensional waveguides 
are prime candidates for such devices. 

Here we introduce a scheme for a quantum single pho- 
ton transistor, a device that can be considered to form 
a cornerstone of quantum optical information process- 
ing. In our approach individual photons propagate in two 
one-dimensional waveguides of low transverse dimension 
and scatter off each other at a localized scattering center 
formed by two two-level systems (qubits) that each cou- 
ple to one of the waveguides, see figure [IJl The qubits 
interact in such a way that no excitations can be ex- 
changed between them and thus ensure that each pho- 
ton remains in its initial waveguide after the scattering 
event. Nonetheless, as we show below, the presence of a 
single photon in one waveguide can completely block or 
enable the propagation of a photon in the other waveg- 
uide. Importantly, our approach considers propagating 
light signals that all have the same carrier frequency and 
is thus inherently scalable as the output signals of one 
transistor can enter as input signals into a consecutive 
transistor, c.f. figure [1J for an illustration. Such scal- 



ability is questionable in previous proposals which are 
based on different technological platforms [3] [4] . More- 
over the device we propose is a passive element that does 
not require any temporal tuning. This implies that the 
arrival time of the photons at the scattering center can 
be completely unknown. 

A technology that is ideally suited for realizing the de- 
vice we envision is provided by itinerant microwave pho- 
tons in superconducting circuits [5j [6]. Here, coherent 
scattering at a superconducting qubit [7H9] and entan- 
glement with a qubit [10] have been demonstrated for 
individual photons that propagate in open transmission 
lines. Moreover precise shaping of single photons pulses 
has been shown [IT] very recently. An implementation of 
our approach in circuit quantum electrodynamics thus 
requires two superconducting qubits that are coupled 
to open transmission lines. We show that the desired 
qubit-qubit interaction can be realized with two trans- 
mon qubits [12] that are coupled via a SQUID which can 
be tuned to ensure that no excitations are exchanged be- 
tween both transmons. Importantly, this coupling is not 
dispersive [13] and thus strong as both transmons have 
the same transition frequency. These rather unique pos- 
sibilities for qubit-qubit interactions offered by supercon- 
ducting circuits are very suitable for our aims. Moreover, 
their robustness with respect to dephasing noise make 
transmons ideal qubits for our device. Yet, alternatively 
one could also use two flux qubits that are coupled via 
an induction loop p~4j [15] . 

To demonstrate the capabilities of the single photon 
transistor we propose, we provide results for the photon 
reflection and transmission spectra and probabilities for 
both transmission lines, which depend on the incoming 
single or two-photon pulses. Furthermore, we determine 
the achievable strength of photon-photon interactions un- 
der realistic experimental conditions, i.e. taking into ac- 
count all dissipative processes in our setup. 

Setup We consider two interacting qubits that each 
couple to a one-dimensional waveguide in which the pho- 
tons propagate. Here we focus on a setup for which we 
can refer to a control and a target photon, where the 
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scribed by the Hamiltonian [20], iJ^ = J^dpp (r^r p — 
I pip) + dppbjjbp, where (Z£) creates a photon in 
subsystem 1 which travels to the right (left) and frj cre- 
ates a photon in subsystem 2. p = v g r] where v g is the 
group velocity and the wave vector n is negative (posi- 
tive) for left (right) going modes. The transmission line of 
subsystem 2 is cut and can be described by only one con- 
tinuum of modes since the incoming and outgoing modes 
of a semi-infinite transmission line can be mapped to an 
infinite waveguide where photons only propagate in one 
direction [T7] . 

The dispersion relation of a transmission line is lin- 
ear and the frequency integration can be extended to 
±00 since we only consider single photon pulses with a 
frequency width that is much smaller than their carrier 
frequency. For these narrow linewidth pulses we thus 
take the photon-qubit coupling to be independent of the 
photon frequency, 



FIG. 1: Sketch of the considered setup, a) Two photons 
traveling in separate waveguides scatter off each other at a 
scattering center formed by two interacting qubits. b) Since 
all photons have the same frequency, multiple devices can 
be concatenated to form a network, c) Level scheme of the 
scattering center Hamiltonian H sys , c.f. equation 



presence of the control photon influences the target pho- 
ton's direction of propagation, while the control photon's 
direction of propagation always changes. A sketch of this 
setup is shown in figure [I] The control photon propagates 
in the waveguide of subsystem 2, c.f. figure [I] which has 
a closed end right where it couples to qubit 2. This ar- 
rangement enhances the absorption of photons by qubit 
2 and hence its inversion as compared to an open waveg- 
uide end, a property that is not necessary but beneficial 
for our aims. The control photon is thus always reflected. 
The target photon in turn propagates in the waveguide of 
subsystem 1. The qubit-qubit interaction is such that no 
excitations are exchanged between the two qubits. This 
implies that a photon in one waveguide can not tunnel 
to the other waveguide and vice versa. Nonetheless one 
control photon in waveguide 2 can completely block re- 
spectively enable the propagation of a target photon in 
waveguide 1 as we show below. The Hamiltonian of the 
two coupled qubits reads, 
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where the <jf are Pauli-operators, uo\ and co>2 the transi- 
tion frequencies of the two qubits and J the strength of 
their mutual interaction. This Hamiltonian can be im- 
plemented with two transmon qubits that are coupled via 
a SQUID, see figure [3] and supplementary material [16], 
or with two inductively coupled flux qubits p~4j [15] . 

Both waveguides respectively transmission lines have 
a continuous spectrum of photonic modes and can be de- 
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Here t\ and T2 are the lifetimes of the two level systems 
associated to their coupling to the transmission lines. A 
further benefit of the semi-infinite waveguide in subsys- 
tem 2 is the y/2 enhancement of the coupling in Hj. Fig- 
ure [TJd shows the level scheme of the two qubits described 
by H sys and the transitions induced by the photons. 

In a realistic system, the qubits will be subject to dissi- 
pation. We thus consider for each qubit a coupling to an 
environment which we model as a continuum of bosonic 
modes with creation operators cj respectively d p and the 
Hamiltonian Hb = dp p c p c p + J^dppdpdp. The 
coupling of the qubits to their environments reads Hbi = 



i) + 
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where 



7 is the relaxation rate of the qubits. The total Hamil- 
tonian that includes the environment, the transmission 
lines, the qubits and all couplings between these con- 
stituents thus reads, 



H = H 7 



Ha 
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(3) 



To investigate the dynamics of single photon pulses in 
this setup, we combine quantum scattering theory [18] 
with the input-output formalism [19] of quantum optics 
as in [20] , see supplementary material [16] for details. 
For the dynamics of the qubit operators, including dissi- 
pative processes, we derive quantum Langevin equations 
[2T] starting from the Hamiltonian ([3]). The solutions of 
these yield the source terms for the input-output rela- 
tions which we then relate to the scattering matrix [16] . 

Photon-Photon Interaction To see the effect of the 
photon-photon interaction most clearly, we first consider 
the situation in which only one target photon but no 
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control photon is present. An incoming target pho- 
ton that travels to the right is described by \Sk t ) = 
J dka t (k)r^(k)\0), where k labels the frequency compo- 
nents. We assume for the target photon a pulse with a 
Lorentzian frequency distribution, a t (k) = {^/7rr t [i(co> t — 
A;)+r t " 1 ]}- 1 . Here r t is the temporal width of the pulse 
and Lj t its carrier frequency. A pulse of this form would 
for example describe a photon that was spontaneously 
emitted into the transmission line from an excited qubit 
as experimentally realized in [22] . We here chose to oper- 
ate the transistor such that the target photon is reflected 
in the absence but unaffected in the presence of the con- 
trol photon and choose u t to be equal to the frequency 
of the transition \g1g2) —> 1 61*72) (\9j) /\ e )j denotes qubit 
j in the ground/excited state), i.e. uj t = <^i + 2 J. The 
reverse mode of operation where the target photon is un- 
affected in the absence and reflected in the presence of the 
control photon can be selected by choosing uu t = uj\ — 2 J 
and works equally well. Without control photon evi- 
dently no photon-photon interaction can take place and 
the output state reads, |* Q ut) = / dp^2 i=rlc A(p)4l )' 
where the transmission amplitude is denoted /3 r (p), the 
reflection amplitude j3i (p) and the amplitude for the tar- 
get photon being lost /3 c (p). These amplitudes relate to 
the initial state via Pi{p) = j dka t (k) Si(k,p), where the 
Si(k,p) are the S-matrix elements for the different pro- 
cesses, see supplementary material [16]. The resulting 
transmission probability for the target photon reads, 
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and the reflection probability pu = r t /(l + T17) (n + r t + 
r t r i7)- We note that Pt + Pr < 1 because the photon 
can also be lost due to qubit relaxation. Importantly, in 
the regime of 7 -1 ^> r t ^> ti, the reflection probability 
for the target photon approaches unity [3]. 

Next we consider the case of the same incident target 
photon but now in the presence of an incoming control 
photon. As the control photon inverts qubit 2, the scat- 
tering center is in the state \g1e2) and the target pho- 
ton can only couple to the transition \g1e2) —> \e1e2), 
see figure ^p. This transition is detuned by 4J from 
the target photon frequency, and thus the transmission 
probability for the target photon approaches unity if J 
is large compared to the linewidths of target pulse and 
qubit 1, J ^> r-f 1 + r t -1 . Our scheme works best if 
the control photon pulse is chosen such that it maxi- 
mally inverts qubit 2. A suitable pulse is thus the time 
reversed version of a pulse resulting from spontaneous 
emission of qubit 2 into the transmission line [23] which 
is often called an inverting pulse [24]. The generation 
of inverting pulses and their release into a transmis- 
sion line was demonstrated recently [11 . For the cut 
transmission line in subsystem 2 an inverting pulse of 
carrier frequency uj c and temporal width r c thus reads 
|* c > = Jdka c (k)b{\0) with a c (k) = {^T c [-i(uj c - 



k) + t c -1 ]} -1 . We note that our results do not change 
if the target photon pulse also has the shape of an in- 
verting pulse so that our scheme is indeed scalable. Due 
to the coupling to vacuum, the qubit 2 is of course never 
completely inverted. The output state reads, |^ ut) = 
Jd Pl d P2 Ei= r ,i,c^j=b,d^M^2)i 1 iPi)jHP2)\0). The 
first index in the amplitudes Pij(pi,P2) refers to the tar- 
get photon, which can be reflected, Z, transmitted, r, or 
lost, c, and the second index refers to the control photon 
which can be reflected, 6, or lost, d. For the probability 
of the target photon being transmitted in the presence of 
a control photon we thus get, 



PTC = dp ± dp 2 |/3 r ,6(Pl,P2)| 2 + |/?r,d(Pl,P2)r 



(5) 



We now quantify the performance of the single pho- 
ton transistor we propose via the difference C s and ratio 
R s between the transmission probabilities for the target 
photon in the presence and absence of a control photon, 



C s Pre ~ Pt and R s Ptc/pt, 



(6) 



where pre and pr are given in equations Q and ^ re- 
spectively. For C s = 1 the setup would describe an ideal 
quantum transistor for single photons. Figure [2] shows 
the achievable transmission contrast, C s , and on-off ratio, 
R s , for a realistic device with given qubit-qubit coupling 
of J = 0.01 uj\ as a function of the qubit decay rate 7. 
As the plot shows, an ideal single photon transistor can 
be realized in the limit of vanishing 7/^1 whereas very 
good performance can already be expected for currently 
realized values of J - 50 MHz and 7/cJi - 10~ 6 [25U27]. 
where a single control photon changes the transmission 
probability of the target photon by a factor 20. 

The performance of the single photon transistor we 
propose depends on the shapes of the target and con- 
trol photon pulses and the parameters of the Hamilto- 
nian ([3|. As expected the best choices for the carrier 
frequencies of the control and target pulses are equal to 
the transition frequencies of ^1^2) —> I #1^2) respectively 
\g1g2) "~ ^ l e i#2)> i- e - oj c = UJ2 + 2 J and uj t = uo\ + 2J. For 
a single transistor oj\ and UJ2 m ay be chosen arbitrarily. 
Yet to enable concatenation of multiple transistors, we 
choose uj\ = UJ2- Moreover the interaction of the target 
photon with qubit 1 should be as high as possible. We 
choose r\ = 20/cji which is compatible with experiments. 
For a control photon which is an inverting pulse, i.e. the 
time reversed version of a photon that would result from 
a spontaneous emission process, the optimal choice for 
its temporal width is obviously r c = T2. There are thus 
two remaining parameters, r t and r c , which we have op- 
timized numerically. The results are shown in figures^ 
and [2)3, where we show the optimal choices of r t and r c 
as a function of 7/^1 . 

Coupled transmons As stated above, the qubit-qubit 
interaction in equation ([T]) can be realized with two 
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FIG. 2: Performance of the quantum single photon transistor, 
a: Maximized contrast C s — pre — Pt as a function of the 
inverse Q-factor of the qubits 7/cJi for uj c = 002 + 2 J, uo t = 
uj\ + 2 J, CJ2 = cji, r\ — 20/cji, J = 0.01 cji and r c = T2. 
b: On-off ratio R s = Ptc/pt for the same parameters, c: 
Optimal choice of rt for a given 7/cJi. d: Optimal choice of 
t c for a given 7/cJi . 



transmon qubits that are coupled via a SQUID. The 
circuit we consider is sketched in figure [3J and de- 
scribed by the Lagrangian [28], C = C\ + £2 + £12, 
where the Lagrangians of the individual qubits read, 
A' = f^j + _ ^') 2 + Ej cos (ipj/(po) and £12 = 

^ (0i - ^2) 2 + ^Jm cos (^™) • Here ^ = fi/(2e) is 
the flux quantum divided by 27r, C and Ej are the ca- 
pacitance and Josephson energy of the individual trans- 
itions which are assumed to be identical. C m and Ej m 
are the capacitance and Josephson energy of the capaci- 
tively shunted coupling SQUID. The C g are the coupling 
capacitances between the transmission lines and the in- 
dividual transmons (also assumed to be identical) and 
the Vi are the fully quantum mechanical quadratures of 
the electric potential of the transmission line fields. All 
Josephson energies of the setup are tunable by threading 
external fluxes <l>j and <£j m through respective SQUID 
loops, c.f. figure [3] We write the corresponding Hamil- 
tonian of the transmons in terms of creation and anni- 
hilation operators at and aj [16 . By tuning Ej and 
Ej m such that ^f+Ej = c+c^+c 7 a ^ interactions 
of the form a\a\ + a\a2 cancel and the leading term of 



(C+C,„+C ff ) a l aia 2 a 2 



the remaining interactions reads — ^?^q C ^ 
which is equivalent to the interaction in equation with 
J — 2{c+c m +c ) Ec- Within the approximations we use 
[T6] the achievable qubit-qubit coupling is J < Eq/ '10. 
In conclusion, we have introduced a scheme for a quan- 
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FIG. 3: Circuit model of the two coupled transmon qubits 
with Josephson energy Ej and shunting capacitance C. Both 
transmons are coupled via a combination of capacitive and 
inductive coupling realized with a SQUID arrangement with 
Josephson energy Ej m and shunting capacitance Cm- The 
transmons are capacitively coupled to transmission lines. 



turn single photon transistor in circuit quantum electro- 
dynamics that is inherently scalable as both photons can 
have the same carrier frequency and pulse shape. The 
performance of the device might be further improved by 
suppressing losses with multiple, regularly spaced qubit 
pairs [29 . Moreover, the complexity of a network built 
with such transistors could be increased further by inte- 
grating directional couplers between them [30] . 
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SUPPLEMENTARY MATERIAL 
Scattering theory and input-output relations 

In this work we calculate the output states of our de- 
vice for given input states by making use of a combination 
of scattering theory and the input-output formalism. In 
this approach, output and input states are connected via 
the scattering matrix (S-matrix), 

|*out> = SI**) (7) 

The S-matrix elements in turn can be written in terms 
of scattering operators [18]. For example for the process 



of transmission of a single photon, the S-matrix element 
reads, 



S r {k,p) = (0\r out (p)rl{k)\0), 



where 

J 



\Hto — \H T to ] ( \ AH T to —iHto 



/° utV ' £ ^Too x ; 
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, (9) 



and H T = dpp (r+r p - lp p ) + dppb\b v is as in 
the main text. In turn, for a two photon process, where 
the target photon is transmitted and the control photon 
is lost in the vacuum, the S-matrix element reads, 

Sr id (k,k' = (0\r ou t 7P d out ^r} nk bl nkf \0) (10) 

For our system it is convenient to introduce even and odd 
modes for transmission line 1, 



1 



ai(p) = -j={r p -l- v ) 



(11) 
(12) 



Due to the form of the photon-qubit coupling in equation 
(2) of the main text, only even modes couple to the qubit 
whereas odd modes completely decouple from the rest of 
the system and describe freely propagating photons. The 
coupling strength to the even modes is enhanced by a 
factor of y/2. 

To calculate the desired scattering matrix elements a 
link to equations of motion, i.e. the Langevin equations 
[T9] that describe the dynamics of the scattering center 
is exploited. To make use of the Langevin equations, it 
is necessary to connect the scattering operators to stan- 
dard input-output operators that fulfill the input-output 
relations, 



(13) 



flout (*) = a>in(t) ~ i\ -a-(t), 



where d[ n (t) and a out (t) are any even mode input /output 
operators, defined by 

am/out (*) = ~j= J oW-^-^ao/iM. (14) 

Here, the operator a /i(o;) plays the role of an initial 
value in the Heisenberg picture, 



a /iM = e 0/1 a(^)e" liit °/ 1 



(15) 



The corresponding input and output scattering operators 
for even modes are defined as in equation (l9|), 



iHTto —iHto 



^in/out M = lim e lHto e- lHTto a(uj)e lHTto e 
= lim e iu;to e iHto a(uj)e- iHto , 



(16) 
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FIG. 4: Circuit model of the pair of transmons with Joseph- 
son energy Ej and shunting capacitance C. Both transmons 
are coupled via a combination of capacitive and inductive cou- 
pling realized with a SQUID arrangement with Josephson en- 
ergy Ejm and shunting capacitance Cm- The transmons are 
capacitively coupled to transmission lines and all Josephson 
energies of the setup are tunable by threading external fluxes 
<E>j and &j m through respective SQUID loops. 



where we have made use of the Baker-Campbell- 



Hausdorff formula. We now solve equation (16) for a{uS) 
and plug the result into equation (14) by making use of 
equation ( [l5| ). Taking the limits to — >• — oo respectively 
t\ — > +oo in equation (14) we find that the scattering 



operator a- m (uj) is the Fourier transform of the Heisen- 
berg operator a- m (t) in the limit to — >• — oo and that the 
the scattering operator a out (oj) is the Fourier transform 
of the Heisenberg operator a out (t) in the limit t\ —> -f oo 



1 / out 



(*) 



= — / 



dua u 



i/out 



-iout 



(17) 



Using equation (17) we express the scattering matrix S 



in terms of input and output operators which are found 
from solutions of the Langevin equations that describe 
the dynamics of the scattering center. 



Coupled Transmons 

The Lagrangian describing the two transmons that are 
mutually coupled via a SQUID and each couple to a 



transmission line reads, 
C 



C ( -2 

2 K 



a 



C m t . 

c„ 



2 
Ej 



■ \ 2 

v 2 f 



(18) 



COS 



Ejm COS 



<P1 



COS 



<£>2 



where cpo = H/(2e) is the flux quantum divided by 27r, 
C and Ej are the capacitance and Josephson energy of 
the individual transmons which are assumed to be iden- 
tical. C m and Ej m are the capacitance and Josephson 
energy of the capacitively shunted coupling SQUID. The 
C g are the coupling capacitances between the transmis- 
sion lines and the individual transmons (also assumed to 
be identical) and the V{ are the fully quantum mechani- 
cal electric potential quadratures of the transmission line 
fields. Hence, the corresponding Hamiltonian reads, 

9 . ~ Cm 



H = 4E C 



8E c n 9jl [ m 
SE c n gy2 

Ej (COS (j. 



c + c g 

Cm, 



-niri2 



(19) 



C- 

Crr, 



■ Cn + C r 



-n 2 



c + c a + c r 



-rix + n 2 



COS 2 ) 



E Jm cos(0i 



where 
E c = 



C+C.+Crr, 



PoCgVi for i 



1,2 and 



The secondary capacitive 



{C+C g ){C+C g +2C rn )' 

coupling between transmission line 1 and transmon 2 and 
vice versa is small compared to the other couplings since 
C m /(C + C g + C m ) <C 1. We thus neglect such cou- 
plings and obtain for the coupling between the transmis- 
sion lines and the transmons, 



H g = SEcirig^ni + n 9j2 n 2 ) 



(20) 



We furthermore separate the Josephson terms into local 
and nonlocal terms, 

Ejm COs(0i - (j) 2 ) = Ejm (COS 01 + COS 2 ) (21) 

s v ' 

local 

+E Jrn [(cos 0i - 1) (cos 02 - 1) + sin 0i sin0 2 ], 



nonlocal 

where we have added an irrelevant constant. This leaves 
us with the local Hamiltonians for each transmon (i — 
1,2), 

Hi = ±E c n\ - Ej cos fc, (22) 
where Ej = Ej + Ej m and the coupling Hamiltonian, 

C m , 



H12 = 8Eq 



C + C Q 



Ejm [(COS 01 



-^-n x n 2 
• C/ rn 

- 1) (COS 02 



(23) 



1) 



sm 0i sm 02j 



7 



We describe the transmons in the approved approxima- 
tion with anharmonic oscillators [12] and introduce rais- 
ing and lowering operators a\ and ai via 



(2Ec 
\Ej 



(a» + aj) 



(24) 
(25) 



Keeping only the leading nonlinear terms we thus find, 

(26) 



Hi 



SE c Ej a\ ai — a\ a\ a^ai 



We are interested in a scenario where tunneling of exci- 
tations from one transmon to the other is strongly sup- 
pressed. The leading tunneling terms are found by ex- 
panding the coupling Hamiltonian H12 to linear order in 
a\ and With a rotating wave approximation we find, 



c + c a + c r 



-nin 2 - Ejmcfixfo 



C -\- C q -f- Cf, 



Ej r> 



(27) 

{a\a\ + a\a2) 



These terms vanish if one chooses the external fluxes that 
control the values of Ej and Ej m such that, 



Ei 



Cn 



Ej + Ej m C + C„ 



C a 



(28) 



The rotating wave approximation applied for deriving 



equation (27) furthermore requires that Ej m < 2Ej/3. 

The desired density-density interaction is contained in 
the nonlocal cosine interaction terms, 



E 



-2E C 



E 



Jm 



Ej + E j7 



-a\aia\a2 (29) 



Here we have applied a rotating wave approximation and 
made use of the fact that our initial states only contain 
one photon per transmission line which allows us to ne- 
glect terms of the form aiaiajaj + H.c. Since the dy- 
namics of our system is, for the initial state we consider, 
restricted to the subspace of at most one excitation per 
transmon respectively transmission line we can write our 
Hamiltonian in terms of Pauli matrices, 



H 



2x2 



He 



H, 



2x2 



(30) 



where H sys is as in equation (1) of the main text with 
j + 5lu : where Suj accumulates all minor 
renormalizations of the transmon frequencies due to the 
nonlinear terms of the cosine potentials and 



J 



2(C + C m + Cg) 



Ea- 



rn 



Since C m / (C + C g + C m ) Cl we find the upper limit, 



J< V0 E C , 



(32) 



for the strength of the qubit-qubit interaction. In turn 
the couplings to the transmission lines read, 



H 



2x2 



which become identical to equation (2) after applying a 
rotating wave approximation. 



